IFIC/10-45 



Neutrino phenomenology and stable dark matter with At 

D.Meloni* 1 , S. Morisi t2 and E. Peinado 12 
1 Institut fur Theoretische Physik und Astrophysik, Universitat Wurzburg, D-97074 Wurzburg. 
2 AHEP Group, Institut de Fisica Corpuscular - C. S. I. C./ Universitat de 
Valencia Edificio Institutos de Paterna, Apt 22085, E-46071 Valencia, Spain 

We present a model based on the A4 non-abelian discrete symmetry leading to a predictive 
five-parameter neutrino mass matrix and providing a stable dark matter candidate. We found an 
interesting correlation among the atmospheric and the reactor angles which predicts #23 ~ 7r/4 
for very small reactor angle and deviation from maximal atmospheric mixing for large #13. Only 
normal neutrino mass spectrum is possible and the effective mass entering the neutrinoless double 
beta decay rate is constrained to be \m ee \ > 4 ■ 10 -4 eV. 



I. INTRODUCTION 

Neutrino oscillation and dark matter (DM) are so far the most important evidences of physics beyond 
the standard model. Many models for neutrino masses and mixings have been studied in literature, some 
of them with the aim of reproducing the tri-bimaximal mixing pattern (TBM) [l| observed, within one 
standard deviation, in the recent neutrino oscillation data [1, E|, namely sin 2 #12 = 1/3, sin 2 6*23 = 1/2 and 
sin 2 = 0. In this context, non-Abelian discrete flavor symmetries have been largely used since, especially 
in models based on A4 [4| and S4 [5[ permutation groups, the TBM limit can be easily reproduced. 

While we have robust evidence for dark matter from rotation curves of spiral galaxies, gravitational lensing, 
WMAP measurament, CMB anisotropy, structure formation, X-ray observations, bullet-clusters [gj], we still 
do not have neither theoretical nor experimental indications about the nature of it; also the mechanism for 
DM stability is still not understood yet. Many models assume ad hoc extra Abelian symmetries 0] in order 
to stabilize the dark matter. Such a symmetry can arise from the spontaneous breaking of a non-Abelian 
continuous symmetry, see for instance or from the breaking of Abelian U(l) group 0- 11 1 as, for example, 
in grand unified frameworks or in supersymmetry [12j |. 

Recently a model where the stability of the DM arises from the spontaneous breakdown of a non-abelian 
discrete flavor symmetry was proposed [13]. Such a model is based on the property that the A4 group is 
spontaneously broken into its subgroup Z2 which is responsably for the stability (for a paper with decaying 
dark matter with discrete non-Abelian symmetry see 0-0]). In |13| ] the same Z2 is also acting in the 



neutrino sector, giving a vanishing reactor angle #13 = and allowing only the inverted hierarchy for the 
neutrino mass spectrum. 

However these properties are model-dependent features and cannot be considered as general results of 
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models where the stability of the dark matter is justified by non-abelian discrete flavor symmetries. In this 
paper we want to provide an explicit example of a model, based on a simple extention of [13| , with a richer 
neutrino phenomenology, predicting a normal mass ordering and #13 7^ 0. 

The paper is organized as follows: in section [TT] we give the field content of our model and derive the 
neutrino mass matrix, whose phenomenological implications are discussed in section IIIII in section IIVI we 
draw our conclusions. 



II. THE MODEL 



In Tab.(|T]) we summarized the model quantum numbers. In contrast to 13 1, the right-handed neutrino N4 
is assigned to 1' instead of 1 and we introduced one more right-handed neutrino N5 assigned to 1". 
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TABLE I: Summary of relevant model quantum numbers 



The operators needed to generate neutrinos masses are the following: 

L = yZLeiNTTih+yZLniNTTih' + y v 3 L T {N T r,) v + 

+y%L T N A H + y^L^N 5 H + M X N T N T + M 2 N 4 N 5 + h.c. 

The scalar and the lepton charged current sectors of our model are the same as in [3] and we refer to that 
paper for all the details. After electroweak symmetry breaking, two of the Higgs doublets contained in the 
A4 triplet r\ do not take vev and we have: 

(H°)=v h ^0, ( V ° 1 )=v v ^0 (r ?2 %>=0. (1) 

As a consequence, the Dirac and Majorana neutrino mass matrices are: 
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Using the See-Saw type-I formula we get the following light neutrino mass matrix 
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where we defined: 
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As discussed in more details in 13j, the group A4 is broken by the vev ~ (1, 0, 0) down to the subgroup 
Z 2 , generated by the diagonal A4 generator S = Diag{l, —1,-1}; the Z 2 symmetry acts on the triples fields 
in the following way: 



N 2 
N 3 



-N 2 , h 2 
-N 3 , h 3 



-h 2 , A 2 
-ha, A 3 



-A 2 . 



(5) 



where h% and Aj are respectively the CP-odd and CP-even component of the Higgs doublet rji for % = 2, 3 
and N 2} 3 are the component of the triplet Nt- This residual symmetry is responsible for the stability of 
the lightest combination of h 2 and h 3 , which is the dark matter candidate. In fact it couples only to heavy 
right-handed neutrinos and not to quarks, supposed to be singlets under A4. Such a scalar dark matter 
candidate is potentially detectable in nuclear recoil experiments 171 Il8|. 



III. NUMERICAL ANALYSIS OF THE FULL MASS MATRIX 



The most general complex and symmetric matrix has twelve independent real parameters which, after 
readsorbing the unphysical phases, reduce to nine. The matrix in eq.© has four complex parameters and 
then five real independent parameters: the moduli of a, b and c and the modulus and phase of the combination 
be + k — de %< ^ d . In terms of them, the matrix in eq. Q is rewritten in the following form: 



m u = ab b 2 de ict,d . (6) 




We can relate three of the five parameters in eq. (JSJ (i.e. a, b and c) with the neutrino masses using the 
invariant quantities of the hermitian matrix M 2 — m v m\, : 

Tr {Ml) = T = a i + 2a 2 (b 2 + c 2 ) + 6 4 + c 4 + 2d 2 = 

2 1 2 1 2 
= m Vl + m V2 + m„ 3 

det(M 2 ) = a 4 (b 2 c 2 -2b cd cos c/) d + d 2 ) 2 = 

222 m\ 
= m Vl m V2 m Va (7) 



- [T 2 - Tr(M 2 M 2 )] = (b 2 c 2 - 2b c d cos <\> d + d 2 ) ■ 

[2a 4 + 2a 2 (b 2 + c 2 ) + b 2 c 2 + 2bcd cos cj) d + d 2 ] 

2/2. 2 \ , 2 2 

= m Vl {m V2 + m V3 ) + m V2 m V3 . 

where the symbols Tr and Det refer to the trace and determinant of a matrix, respectively. The expressions of 
a, b and c in terms of neutrino masses, derived from the inversion of the previous relations, produce lengthly 
formulae that we do not present here but included in our numerical simulations. As independent parameters, 
we decide to take m Vl , the two square mass differences Am^, and Am 2 tm and the other two parameters of the 
neutrino mass matrix, i.e. d and cf>d, to study the correlations among the neutrino observables predicted by 
the model, we perfom a parameter space scan extracting randomly the previous parameters in the following 
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intervals: 

de[-l,l], (f> d e [-7r,7r), m Vl £ [10- 5 , 1] eV, 
7.1 x 1CT 5 eV 2 < Am 2 sol < 8.3 x 1(T 5 eV 2 , 

0.027 < r = ,£ m °° l < 0.040. 

Then, we require that the obtained mixing angles are within their current 3er confidence level Q. We 
carefully checked that \d\ > 1 is not compatible with data. We made the calculation for both inverted 
and normal neutrino mass hierarchies finding that only the latter is allowed in our framework (we have 
analitically discussed this point in the appendices: Appendix IA II for the real case, Appendix IA 21 for the 
the fj, — t invariant case and Appendix IA 31 for the exact TBM case). Among all possible correlations, we 
found that the (#13 — #23) one is quite interesting, as shown in fig- (HJ - As we can see, for very small #13 the 




atmospheric angle is fixed to its maximal value, as a consequence of the fact that the model has the TBM 
limit built-in (see next section). Increasing #13 two different branches above and below 623 = 7r/4 develop 
and for #13 J> 1° maximal 2 — 3 mixing is strongly excluded. It is interesting to observe that this result can 
be easily verified at the forthcoming neutrino experiments 21] even with a reduced sensitivity to #13 since 
the largest deviation from 823 — 7r/4 is obtained for #13 close to its current upper bound. We checked that 
the previous correlation is obtained with a, b, c and d all at the same order of mag nitude 0(1O -2 ), which 
ensures that no fine-tuning is required among the matrix elements of eq. (J3J) , and with the phase <f> uniformly 
distributed in the [— n, n) interval. We also verified that the model does not produce any other relevant 
correlations among the mixing parameters. 

The other interesting neutrino observable that we want to discuss is the effective mass \m ee \ entering 
the neutrinoless double beta decay rate. Since our model is only compatible with the normal hierarchy, we 
expect it to be quite small for small m Vl . The result of our simulation can be found in fig. ([2]), where we 
can see that a lower bound (\m Vl |, |m ee |) *S (2 • 10 -3 , 4 ■ 10~ 4 ) eV can be set. The existence of such a lower 



5 



bound on \m ee \ can be justified in the following way: since we are working in the basis where the charged 
leptons are diagonal, a vanishing \m ee \ would imply a vanishing parameter a; in that case, the mass matrix 
in eq. (J6j) would have a zero mass eigenvalue associated to the eigenvector (1, 0, 0) T , which is not compatible 
with the neutrino experimental data. 




FIG. 2: Values of the effective mass \m ee \ (in eV) allowed in our model. The two dashed horizontal lines represent 
the experimental sensitivities of some of the forthcoming experiments while the dashed vertical line is the upper limit 
from tritium f3 -decay experiment. For references to experiments see jjl - l^} /. 



IV. CONCLUSION 



In this paper we have studied a model explaining the stability of dark and giving an interesting neutrino 
phenomenology based on the A4 flavor symmetry. The model is an extension of the standard model and 
contains four Higgs doublets and five heavy right-handed neutrinos. Light neutrino masses are generated 
with the type-I seesaw and the resulting Majorana mass matrix has only five free parameters. We made 
a numerical scan of the allowed parameter space and found a correlation between the atmospheric and the 
reactor angles. In particular, for #13 J> 1° maximal atmospheric angle is strongly disfavoured whereas for 
small reactor angle the #23 is close to maximal. The model gives only normal neutrino mass hierarchy and 
predicts a lower bound for the neutrinoless double beta decay \m ee \ ^ 4 • 1CT 4 eV. 
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Appendix A: Mass hierarchy 



1. Neutrino mass hierarchies for the real mass matrix 



In this section we show that our model allows only the normal hierarchy in the simplified case the mass 
matrix in eq.© is real. Since the absolute neutrino mass scale is unknown, we can restrict ourselves to study 
the situation where the lightest neutrino mass eigenvalue is (almost) vanishing. For our discussion it is then 
enough to consider the determinant of such a matrix, which results: 

det(m„) = -a 2 {-bc + d) 2 . (Al) 

For an almost vanishing determinant, we can have two options, namely a ~ and d ~ be. In the first case, 
the neutrino mass matrix allows a diagonalizing matrix of the form: 

/ 1 \ 

m v ~ x z , (A2) 
\0 y f J 

whose first line implies a vanishing solar mixing angle. The phenomelogical interesting case corresponds then 
to a but, even in this case, the null eigenvalue will always be associated with an eigenvector different 
from (0, 1, — 1) T , so it cannot be m Ui anyway. This excludes the inverted hierarchy as a viable neutrino mass 
spectrum. The second choice d = be always produces two vanishing neutrino masses; if one of them would 
be m V3 then the atmospheric mass difference would be vanishing; on the other hand, if Wj VX — 0, then one 
can have Am^i = but Am^ 7^ 0. Then the normal mass scheme provides the only framework where to 
account for such results. 



2. Neutrino mass hierarchies for the vanishing reactor angle 



We work here in the limit of vanishing #13, which is a good approximation as given by the experimental 
data. In our case this limit implies b — c m eq. This matrix has an eigenvector (0, —1, 1) with eigenvalue 

m„ 3 = (c 2 - de***). (A3) 

The absolute value of this mass squared is 

\m„ 5 \ 2 = c 4 - 2c 2 d cos (j) d + d 2 . (A4) 
From the invariant equations ([7|. the determinant of the squared mass matrix is 

det(M 2 ) = a 4 (c 4 - 2c 2 dcos d + d 2 f = \m Vl \ 2 |m„ 2 \ 2 \m V3 \ 2 . (A5) 
From eqs. (IA4|) and (|A5[) we have the relation 
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we also know that |to&!,| — a 2 and therefore 



\ m 12 / A_2 



|mbfc| = Ki1 krF = Ki1 I 4 + 4 + 4 r# + 44 bfcn — cosa J (A7) 

where a is the Majorana phase and c© = cos 6*12 and Sq = sin #12. From the previous relation we get: 



feP = 4 + 4 + 4gg| + c|4 ^g^g cos. ) ( As ) 



which, for |m v x| >> V Araf 2 , can be approximated by: 

~ ( c © + s © + c © s © C0Sa ) = 1 _ C © S Q (2 _ C0Sa) < 1 (A9) 

That means that |?n V 2 1/1^1/31 < 1 and implies a normal hierarchy. In the limit m v \ ~ AmJ 2 we cannot 
have the inverse hierarchy because the minimal value for m„i is m™ n = yAm^ > y Amf 2 - 

Notice that the relation (|A8I) can be expressed as an implicit function of m^i, Amf 2) Am^, #12 and the 
Majorana phase a. Using the experimental information on the mass differences and the solar angle we can 
numerically evaluate the minimum allowed value for m v \ . 



3. TBM limit 



In this appendix we briefly discuss the TBM limit of the mass matrix in eq.([3]). This can be achieved 
imposing the following relations among the parameters: 

b = c 

(A10) 

k = a 2 + ab- 2b 2 

(the first relation is enough to get a /1 — r invariant mass matrix) . Then the masses depend on two complex 
parameters that can be ordered according to: 

mi = a(a — b) 
777,2 = a(a + 2b) 

m 3 = -a 2 - ab + 2b 2 . (All) 

It is easier to study the phenomenology redifining 

a = \a\e 1 ^ 
b = l^le^ 6 

H = t {4> a - fo) = M 

\a\ 

so that the solar and atmospheric mass differences are: 
A77721 = 3|a| 4 t(2cosA0 + t) 



(A12) 



Amjj = 4|a| 4 t(cosA0 + /j- 2/jcos^ A0-r cosAt/f + r). 
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The model is only compatible with a normal hierarchy spectrum because the simultaneous requirements 
Amli > and Am| x < gives 



t 



< cosA0 < 1 



cosA0 < — t 



(A13) 



which are obviously incompatible. Moreover, it easy to check that the conditions \mi\ > and 1 777.3 1 ^0.5 
eV imply: 



0.07 < t < 5 . 



(A14) 



For |77i ee | we found a lower bound |777 ee | > 6 • 10 4 eV; other neutrino mass matrix models with two complex 
parameters predicts different lower limits in the TBM limit, for instance |m ee | > 7 • 10~ 3 eV in [27J and 
\m ee \ = in [28 1 . 
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